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Abstract 

In the minimal supersymmetric model (MSSM) all Higgs self-coupling param- 
eters are related to gauge couplings at tree-level. Leading-logarithmic radiative 
corrections to these quantities can be summed using renormalization group tech- 
niques. By this procedure we obtain complete leading-log radiative corrections to 
the Higgs masses, the CP-even Higgs mixing angle, and trilinear Higgs couplings. 
Additional corrections due to squark mixing can be explicitly incorporated into 
this formalism. These results incorporate nearly all potentially large corrections. 
Mass shifts to the neutral CP-even Higgs bosons grow with the fourth power of 
the top-quark mass and can be significant. The phenomenological consequences of 
these results are examined. 
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1. Introduction 



Supersymmetry is one of the most promising theoretical ideas that attempts to 
explain the origin of the scale of electroweak interactions. The minimal supersym- 
metric extension of the Standard Model (MSSM) is the most economical among 
models of this type [1], and deserves close examination as a candidate for a model 
of physics beyond the Standard Model (SM). In the MSSM, one simply adds a 
super symmetric partner to every quark, lepton and gauge boson. In addition, the 
MSSM must possess two Higgs doublets in order to give masses to up and down 
type fermions in a manner consistent with supersymmetry (and to avoid gauge 
anomalies introduced by the fermionic superpartners of the Higgs bosons). The 
Higgs sector of the MSSM is greatly constrained by supersymmetry [2]. All quar- 
tic Higgs coupling constants are related to electroweak gauge coupling constants, 
which imposes various restrictions on the tree-level Higgs masses and couplings. In 
particular, all tree-level Higgs parameters can be expressed in terms of one physical 
Higgs mass and the ratio of vacuum expectation values, tan/3 = vijv\. 

Any realistic supersymmetric model must incorporate supersymmetry break- 
ing in the low-energy theory. This breaking is parametrized by adding soft su- 
persymmetry breaking mass terms for the squarks, sleptons and gauginos, and 
trilinear Higgs-squark-squark and Higgs-slepton-slepton interactions which are pro- 
portional to the so-called A parameters [3]. Due to these supersymmetry breaking 
terms, one finds that the tree-level relations among Higgs masses and couplings 
acquire radiative corrections. It has been shown that these corrections can indeed 
become very substantial if the top quark mass is much larger than mz [4-20]. 
For example, the tree-level bound, m^o < mz receives a radiative correction of 
order g\m\ jm\ ln(M~ /ml) which raises the upper limit of m^o by as much as 
20 (50) GeV for a top-quark mass of mj = 150 (200) GeV [5]. The logarithmic 
enhancement factor is a remnant of the cancellation of divergences generated by 
virtual particles (in this case the top quark and its supersymmetric partners). Sim- 
ilar logarithmic corrections also arise when contributions from other sectors of the 
theory are incorporated. 

While the exact one-loop radiative corrections to mass sum rules and Higgs 
mass bounds can be obtained in a straightforward manner, the radiative correc- 
tions to individual CP-even Higgs masses and Higgs interactions are far more 
complex. For instance, computations of the latter type require a careful definition 
of the parameter tan f3. However, significant simplification can be achieved if we 
include only the leading logarithmic radiative corrections. In this approximation, 
the definitions of tan/5 and the CP-even Higgs mixing angle a are unambiguous 
and can easily be related to physical observables. The goal of this paper is to 
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construct a low-energy effective Lagrangian from which one can directly obtain the 
leading contributions to the radiatively corrected Higgs masses and couplings. 

One possible procedure for achieving this goal is to compute the terms of the 
full one-loop effective action that depend on the scalar Higgs fields. This requires 
the computation of the effective potential, V e g , the coefficient of the scalar kinetic 
energy term, Z e g, and higher derivative terms. Let us assume that the supersym- 
metry breaking scale (Msusy) is somewhat higher than the electroweak breaking 
scale. Here, we use Msusy to denote the mass scale that characterizes supersym- 
metry breaking. (For now, we ignore the possibility of multiple supersymmetric 
thresholds, which will be addressed later in this paper.) If one expands the ef- 
fective action about the Higgs vacuum expectation value and discards all terms 
of order ^|/Mg USY , then ^ suffices to keep only those terms of dimension 4 or 
less. Note that at this stage, it is not strictly correct to use the effective potential 
to compute the one-loop Higgs masses and couplings. One must first rescale the 
scalar fields (i.e., wave function renormalization) in order that the scalar kinetic 
energy terms are canonical. The end result is a scalar potential that is polynomial 
in the scalar fields (with terms of dimension 4 or less), whose coefficients reflect 
the one-loop radiative corrections. The leading one-loop correction terms will de- 
pend logarithmically on Msusy > which suggests that one can use renormalization 
group methods to explicitly identify the leading logarithmic terms [7,8,13,19]. In 
this paper, we shall employ the renormalization group method for two reasons: 
(i) simplicity, and (ii) the integration of the one-loop renormalization group equa- 
tions (RGEs) effectively sums the leading logarithmic contributions to all orders 
in perturbation theory. This method can be extended to include the effects of 
supersymmetric thresholds. In this paper, we derive RGEs that incorporate terms 
which are logarithmic in the ratio of threshold masses to one-loop. However, to 
sum such effects to all orders significantly complicates the analysis and is beyond 
the scope of this paper. 

In addition to neatly summarizing the most important radiative corrections, 
the direct identification of the terms logarithmic in Msusy P rov ides an important 
check for a more precise (and hence more complicated) explicit one-loop compu- 
tation. Moreover, the full RGE analysis yields leading log terms to all orders in 
perturbation theory, and hence provides some information on the contributions 
that lie beyond the one-loop approximation [10]. This provides a check on the 
reliability of the one-loop results. However, it is important to realize that the 



* If Msusy is roughly equal to the scale of electroweak symmetry breaking, then the size of the 
radiative corrections discussed in this paper will be rather small. In this case, one will need 
exact one-loop computations to determine reliably the effects of the radiative corrections. 
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renormalization group technique may not detect all significant terms of the one- 
loop radiative corrections. For example, if the A-terms that control squark mixing 
are large, then radiative correction terms of order A 2 /M| USY may be enhanced by 
powers of the top-quark mass [9,12,17,18,20]. Such terms can then compete with 
(and may be more important than) the leading log terms identified above. Thus, 
an important goal of this paper is to demonstrate how to include such terms within 
the renormalization group approach in a consistent manner. The end result of our 
work is a simple and powerful technique that includes in a transparent fashion the 
most important radiative corrections to the MSSM Higgs sector. 

In section 2, we first examine the general non-supersymmetric two Higgs dou- 
blet model at tree-level. The most general two-Higgs-doublet model potential 
depends on three mass parameters and seven dimensionless couplings. The Higgs 
mass matrices and three-point Higgs vertices are obtained in terms of these param- 
eters. In section 3, we specialize to the MSSM. Radiative corrections to the Higgs 
sector parameters are obtained by renormalization group evolution. Supersymme- 
try implies definite relations among the Higgs self-couplings in terms of the gauge 
couplings. These relations are used as boundary conditions for the renormalization 
group equations at a scale MgusY which characterizes the scale of supersymmetric 
particle masses. The Higgs couplings are then evolved according to renormalization 
group equations (RGEs) and the radiatively corrected Higgs masses and couplings 
are computed in terms of parameters at the electroweak scale. In the evolution of 
couplings, we first assume that the effective low-energy theory at the electroweak 
scale contains two light Higgs doublets. Another possibility is to assume that the 
second Higgs doublet is significantly heavier than the electroweak scale. In this 
case the low-energy effective theory is identical to the SM with one physical Higgs 
boson. This case is discussed section 4. In our analysis, an important parameter is 
tan {3, the ratio of vacuum expectation values. The relation of tan (5 to physically 
measurable observables is elucidated in section 5. 

The analysis described above is equivalent to summing the leading log radiative 
corrections to all orders in perturbation theory. This approximation summarizes 
almost all potentially large radiative corrections. However, if squark mixing effects 
are substantial, an additional set of non- logarithmic corrections may be impor- 
tant. We incorporate these terms into our analysis in section 6. Numerical results 
and their phenomenological implications are given in section 7. The importance of 
using RGE-improved results for the Higgs masses is subject of section 8. Final con- 
clusions are given in section 9. Details of the RGE analysis and the incorporation 
of important non-leading log terms are relegated to four appendices. 

Some results of this paper were first described in ref. [13] and used in the 
phenomenological analyses presented in refs. [21-23]. Additional details can be 
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found in ref. [24]. 



2. The General Two-Higgs-Doublet Model 

We begin with a brief review of the general (non-supersymmetric) two-Higgs 
doublet extension of the Standard Model [2]. Let $i and $2 denote two complex 
Y — 1, SU(2)i doublet scalar fields. We introduce the notation 



Hi \ 



The most general gauge invariant scalar potential is given by 

V = mf^l^i + m22$2 $ 2 - [mi2$l$2 + h.c] 

+ iAi^i) 2 + ±A 2 ($j,$ 2 ) 2 + M^i){^ 2 ) + A4($i$ 2 )($2^i) (2.2) 

+ {iA 5 ($i$2) 2 + [A 6 ($l$i) + A 7 ($5$ 2 )]$l$2 + h.c.} . 

In most discussions of two-Higgs-doublet models, the terms proportional to A6 and 
A7 are absent. This can be achieved by imposing a discrete symmetry $1 — > — $1 
on the model. Such a symmetry would also require toi 2 = unless we allow a soft 
violation of this discrete symmetry by dimension- two terms. For the moment, we 
will refrain from setting any of the coefficients in eq. (2.2) to zero. In principle, 
m 2 2 , A5, A6 and A7 can be complex. In this paper, we shall ignore the possibility 
of CP-violating effects in the Higgs sector by choosing all coefficients in eq. (2.2) 
to be real. The fields will develop non-zero vacuum expectation values (VEVs) if 
the mass matrix mf- has at least one negative eigenvalue. Imposing CP invariance 
and U(1) EM gauge symmetry, the minimum of the potential is 



W = -f=\ , (2-3) 





where the Vi can be chosen to be real. The VEVs have been normalized so that 
m W = j92( v l + v 2)- ft is convenient to introduce the following notation: 



.2 _ 2 , 2 



v 



vf + vi, t p = tan = v 2 /vi . (2.4) 



The gauge symmetry will then be broken spontaneously. As a result three of the 
eight degrees of freedom of the original Higgs doublets are eaten by the W ± and Z. 



* This latter requirement is sufficient to guarantee the absence of Higgs-mediated tree-level 
flavor changing neutral currents. 
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The remaining five physical Higgs particles are: two CP-even scalars (H° and h°, 
with m/jo < mjjo), one CP-odd scalar (A ) and a charged Higgs pair (H^. The 
mass parameters ran and r«22 can be eliminated by imposing the minimization 
conditions 

m 2 u - tpm\ 2 + \v 2 c\ (\ x + 3A 6 t/3 + \t% + X 7 tfj = , 
m\ 2 - t- p l m\ 2 + \v 2 s\ (\ 2 + SAt^ 1 + X S t p 2 + A 6 t^ 3 ) = , 

where we have introduced the following abbreviations: sp = sin/5 and cp = cos (3 
and 

A 3 = A 3 + A 4 + A 5 . (2.6) 

It then follows that the mass matrices of the CP-even, CP-odd and the charged 
scalars are given by 

M 2 A0 = [m\ 2 - \v 2 {2\5spcp + \ 6 c} + X 7 s})] ^ ^ ^ 



M H ±= [m 12 - \v (A 4 s /3 c /3 + X^spcp + \ 6 c p + \ 7 Sp)] 



r 1 



M 2 H o = m\ 2 | i + _ 1 



-l / 



P ' (2.7) 



+ \v 2 spcp 



2Ait /3 1 + 3A 6 - \ 7 t 2 p 2A 3 + 3(A 6 t /3 1 + \ 7 tp) \ 
2A3 + 3(A6t^ 1 + \7tf3) 2A2t/3 + 3A7 — \§t^ 2 J 



The first two mass matrices possess a zero eigenvalue corresponding to the Gold- 
stone bosons (G°,G' ± ). The masses of the physical Higgs particles are given by 
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m 2 A0 = tr(M 2 A0 ) = — ^ - ^ 2 (2A 5 + A 6 t^ + \ 7 t p ) , 

"m 2 ' (2 ' 8) 
m 2 H± = tr(M 2 H ±) = ^ - \v 2 {\ A + A 5 + Ae^ 1 + \ 7 tp) . 



If we now substitute the remaining mass parameter m\ 2 in favor of m 2 A0 we find the 
following expressions for the charged Higgs mass and the neutral CP-even Higgs 
boson mass matrix 

m 2 H ±=m 2 Ai) + ±v 2 (\ 5 -\ 4: ), (2.9) 
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M 2 H o = m\ 



s 2 p -spcp 

SpCp C 2 



\ic 2 p + 2\ 6 spcp + A 5 s^ (A 3 + A 4 )s /3 c /3 + A 6 c| + A 7 s^ \ 
(A3 + M)sj3Cp + X 6 c 2 + A 7 s| A 2 s| + 2A 7 s /3 c /3 + A 5 c| / 



• (2.10) 



The CP-even Higgs mass eigenvalues are given by 



m 2 H o 



h0 = i|tr7W 2 ±y / [tr7W 2 ] 2 -4detM 2 |, (2.11) 



where M 2 = M 2 HQ and the mixing angle a is obtained from 



sin 2a = = 



V^tr^]^4detA^ 
.M '7, .Vlr' 



cos 2a- ' 11 22 



V[trA^ 2 ] 2 -4detA^ 2 



The phenomenology of the two-Higgs doublet model depends in detail on the 
various couplings of the Higgs bosons to gauge bosons, Higgs bosons and fermions. 
The Higgs couplings to gauge bosons follow from gauge invariance and are thus 
model independent. Most of these couplings are proportional to either sin(/3 — a) 
or cos(/3 — a). In contrast, the Higgs couplings to fermions are model dependent, 
although their form is often constrained by discrete symmetries that are imposed 
in order to avoid tree-level flavor changing neutral currents mediated by Higgs 
exchange [25] . An example of a model that respects this constraint is one in which 
one Higgs doublet (before symmetry breaking) couples exclusively to down-type 
fermions and the other Higgs doublet couples exclusively to up-type fermions. This 
is the pattern of couplings found in the minimal supersymmetric model (MSSM). 
Detailed Feynman rules can be found in ref. [2]. Finally, the 3-point and 4-point 
Higgs self-couplings depend on the two-Higgs-doublet potential [eq. (2.2)]. The 
Feynman rules for the most important trilinear Higgs vertices are listed below: 
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9h°A°A° = ~ [ A l s ^ c P s a ~ ^CpSpCa ~ ^(spCa - c|s a ) + 2\^Sp_ a 

^2 

- \sp{c(3S a+ p + s a c 2 p) - hcp(c a C2/3 - sps a+ p)] , 

2iTn\\r r 2 \ 2 T / 3 3 \ r> \ 

g H o A o A o = [\\SpCpCa + \2CpSpSa + \z{SpS a + CpC a ) - 2A 5 Cp^ a 

$2 

- \sp{cf3C a+ (3 + C a C2/?) + A 7 C /3 (s / 3C Q+/ 3 + S a C 2(9 )] , 

— 6m^ r 2 x 2 T / 3 3 2 \ (2-13) 

9H°h°h° = [MSaCpCa + A 2 C Q S /3 S Q + A 3 (s Q Sa + C a C(3 - gC^_ Q ) 

^?2 

- KSa(cpC 2a + C Q C a+/3 ) + A 7 C Q (s^a + SqC q+/3 )] , 

9h°H+H- = 9h°A°A ~l As " M)Cf3-a , 

92 

2m w , , 

9h°H+H- = 9h°A°A" —\ A5 ~ A 4j s /3-a ■ 

92 

(In our notation, g\ = and g> 2 = <?.) It is interesting to note that couplings of the 
charged Higgs bosons satisfy relations analogous to that of m H ± given in eq. (2.9). 

At present, some experimental constraints on the parameters of the two-Higgs 
doublet model have been obtained at LEP. Here we briefly summarize the results 
for the Higgs search as compiled by the Particle Data Group [26] . For the charged 
Higgs boson, mjj± > 41.7 GeV. This is the most model independent bound and 
assumes only that the decays dominantly into r + u T , cs and do. The LEP 
limits on the masses of h° and A are obtained by searching simultaneously for 
Z -> h°ff and Z -> h°A° [27,28]. The ZZh° and Zh°A° couplings which govern 
these two decay rates are proportional to sin(/3 — a) and cos(/3 — a), respectively. 
Thus, one can use the LEP data to deduce limits on m^o and m^o as a function of 
sin(/3 — a) [28]. Stronger limits can be obtained in the MSSM where sin(/3 — a) is 
fixed by other model parameters. The present limits as summarized by the Particle 
Data Group [26] are m^o > 29 GeV and m A o > 12 GeV based on supersymmetric 
tree-level relations among Higgs parameters, but with no assumption for the value 
of tan/3. If leading log radiative corrections are incorporated and tan/3 > 1 is 
assumed, then recent results of the ALEPH Collaboration yield m^o > 41 GeV 
and m A o > 20 GeV (at 95% CL). However, as was shown in ref. [29] (and will 
be discussed briefly in section 7), the limit on m^o may be substantially weaker if 
large squark mixing is permitted. 

The experimental information on the parameter tan f3 is quite meager. For def- 
initeness, we shall assume that the Higgs-fermion couplings are specified as in the 
MSSM. In the Standard Model, the Higgs coupling to top quarks is proportional 
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to g2mt/2mw, and is therefore the strongest of all Higgs-fermion couplings. For 
tan/3 < 1, the Higgs couplings to top-quarks in the two-Higgs-doublet model dis- 
cussed above are further enhanced by a factor of 1/ tan/?. As a result, some weak 
experimental limits on tan/3 exist based on the non-observation of virtual effects 
involving the H~tb coupling. Clearly, such limits depend both on m#± and tan (3. 
For example, for mjj± ~ mw, limits from the analysis of B°-B° mixing imply 
that tan/3 ^ 0.5 [30]. No comparable limits exist based on top-quark couplings to 
neutral Higgs bosons. 

Theoretical constraints on tan (3 are also useful. If tan (3 becomes too small, 
then the Higgs coupling to top quarks becomes strong. In this case, the tree- 
unitarity of processes involving the Higgs-top quark Yukawa coupling is violated. 
Perhaps this should not be regarded as a theoretical defect, although it does render 
any perturbative analysis unreliable. A rough lower bound advocated by ref. [30], 
tan/3 ^ rn^/600 GeV, corresponds to a Higgs-top quark coupling in the perturba- 
tive region. A similar argument involving the Higgs-bottom quark coupling would 
yield tan/3 ^ 120. A more solid theoretical constraint is based on the requirement 
that Higgs-fermion couplings remain finite when running from the electroweak scale 
to some large energy scale A [31,32]. Beyond A, one assumes that new physics en- 
ters. The limits on tan f3 depend on mt and the choice of the high energy scale 
A [13,31,32]. For example, if there is no new physics (other than perhaps minimal 
super symmetry) below the grand unification scale of 10 16 GeV, then based on the 
CDF limit [33] of m t > 91 GeV, one would conclude that 0.5 < tan/3 < 50. Fi- 
nally, it is interesting to note that these limits on tan f3 are not very different from 
those that emerge from models of low-energy supersymmetry based on supergravity 
which strongly favor tan/3 > 1 [34]. 



3. The Radiatively Corrected Higgs Sector of the MSSM 



In a general two-Higgs-doublet model none of the relations derived in section 2 
are very predictive due to the large number of unknown parameters. However, in 
the MSSM, supersymmetry (SUSY) implies constraints among these parameters, 
thereby leading to numerous predictions for Higgs masses and coupling constants 
in terms of a few basic model parameters. 

Consider first the case of unbroken SUSY. Here all the Higgs self-coupling 
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constants are related to the gauge coupling constants 

Ai = A 2 = i(^ 2 2 + ^), 

A4 = -2^2 > 

A5 = A6 = A7 = . 

As a result the Higgs sector of the MSSM is completely determined by two new 
measurable quantities, which can be conveniently chosen to be m^o and tan f3 [2] . 

However, the parameters of any theory will in general depend on the energy 
scale (y/s) at which they are evaluated. This dependence is described by the 
renormalization group equations (RGEs) 

dpi/dt = Pi(pi,p2, ••) , where t = ln(s) . (3.2) 

Here the parameters pi stand for the Yukawa couplings hf (f — u,d,£), the gauge 
couplings of the SU(3) C x SU(2)l x U(1)y gauge group g{ (i = 1,2,3), the Higgs 
self-coupling constants Xj (j = 1,...,7), and the mass parameters of the Higgs 
bosons mjj = 1, 2). Since eq. (3.1) is valid at an arbitrary energy scale y/s we 
find analogous relations for the corresponding /3-functions (/3 Pi = (3/) by taking the 
derivatives with respect to t 

P\ 3 = ?lPgt-P 9 l}, (3-3) 
(3\ 4 = ■ 

Note that eq. (3.3) is valid only if the theory is supersymmetric at the scale y/s. 
However, if SUSY-breaking terms are included, the mass degeneracy between the 
particles and their supersymmetric partners is violated, and the masses of the 
superparticles can become heavy. For simplicity we will typically assume that 
the supersymmetric particle masses are roughly of the same order. That is, the 
scale of SUSY-breaking is characterized by one single parameter, MgusY- Then 
the /^-functions at an intermediate scale y/s (for mz < y/s < Msusy) wm 110 
longer satisfy eq. (3.3) and the gauge coupling constants and the self-coupling 
constants will evolve differently. In the case of multiple SUSY-breaking scales, as 
y/s decreases one would have to modify the /9-functions every time a multiplet of 
supersymmetric particles decouples. In this section we shall assume that the mass 
parameters of the Higgs potential are of the order of mz- This guarantees that 
both Higgs doublets are present in the low effective theory which is equivalent to 
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(3.1) 



the (non-supersymmetric) SM with an extended two-doublet Higgs sector. The 
required /3-functions are presented in Appendix A. 

We begin by running the gauge coupling constants from the electroweak scale 
(where they are measured) up to Msusy- At this scale, the SUSY boundary con- 
ditions given by eq. (3.1) can be imposed, and we obtain the values of the Higgs 
self coupling constants Aj at MgusY- Next we determine Ai(M wea ^) by integrat- 
ing the corresponding RGEs from Msusy down to M wca k. Finally, using Aj in 
eqs. (2.9)-(2.13), we find the RGE-improved Higgs masses and trilinear interac- 
tions.^ RGE-improved Higgs masses based on the two-Higgs doublet RGEs have 
also recently appeared in ref. [19]. Analytic approximations to the trilinear Higgs 
interactions can also be found in the literature. In ref. [12], the h°A°A° coupling 
was computed using the effective potential method (in which only terms explicitly 
proportional to mf were kept). All other three-Higgs interactions were obtained 
using the same approximation scheme in ref. [20]. In ref. [23], the h°A°A° coupling 
was obtained by the method outlined above, in which the RGE-derived expres- 
sions for \i(M wea k) were inserted into eq. (2.13), and is contrasted with the results 
of the effective potential technique. Since this work was completed, a number of 
more complete computations based on one-loop vertex corrections to the Higgs 
self-couplings have appeared. The correction to the hhh coupling can be found in 
refs. [35] and [36]; see also ref. [37]. 

By using the running parameters of the theory evaluated at the electroweak 
scale (MweajJ, one has incorporated the leading logarithmic radiative corrections 
to the Higgs parameters, summed to all orders in perturbation theory. The RGEs 
can be solved by numerical analysis using the computer. But it is instructive to 
solve the RGEs iteratively. To first approximation we can take the right hand side 
of eq. (3.2) to be independent of ln(s). That is, we evaluate the [3{ by imposing 
tree-level relations among the parameters pi [i.e., eq. (3.1)] and evaluating the 
results at the scale y/s = M wea k. Then, integration of the RGEs is trivial, and we 
obtain 

Pi(s 1 )= P i(s 2 )-fr\n(j?J , (3.4) 
under the assumption that the particle content of the effective low-energy theory 



* In refs. [7-10], it was assumed that to^o ~ .Msusy, in which case the effective low-energy 
theory consists of the non-supersymmetric SM with one physical Higgs boson. Our approach 
extends the results of these authors by allowing for the possibility that all five physical Higgs 
states (h°,H°,A° and i? ± ) may have masses substantially below Msusy- 

f In the leading logarithmic approximation, the masses obtained in this manner are physical 
masses (corresponding to the pole of the Higgs propagator). Similarly, all one- loop def- 
initions of the three and four-point couplings differ only in their non-leading logarithmic 
terms. 
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does not change within the range of integration. 

The lower limit of integration is the electroweak scale, M wea k. One possible 
choice for this scale is M wea k = mz (which is the appropriate scale for diagrams 
involving neutral gauge and Higgs bosons inside the loops). A second equally rea- 
sonable choice would be M weg ^ = m t (which is the appropriate scale for diagrams 
involving the top-quark). This is a somewhat arbitrary decision, since a different 
choice would yield results that differ formally from eq. (3.4) by a non-leading log- 
arithmic term (i.e., a term that does not grow as ln(s2) where ^/s2 is the large 
scale) . On the basis of more precise one- loop calculations, we have adopted the fol- 
lowing strategy. We integrate the RGEs from MgusY down to m t . At that point, 
we formally integrate out the top-quark from the low-energy theory, and finally 
integrate the appropriate RGEs of the new low-energy effective theory down to 
mz- The first order solution of eq. (3.2) becomes 

Pi {ml) = p 4 (M| USY ) - A In (^%) - 01 In (4) , (3-5) 

where the are the /^-functions of the SM with two Higgs doublets presented in 
Appendix A. The /3f are the /3-functions in the same model with the top-quark 
decoupled as we now explain. 

We obtain from /3\. by setting the top-quark Yukawa coupling to zero. 
For /3®2, the situation is more subtle. The reason is that below the top quark- 
threshold, the SU(2)l x U(l)y gauge symmetry is broken. As a result, the cou- 
pling constants for vertices involving the gauge bosons are no longer constrained 
by the SU(2)l x U(1)y gauge symmetry and can evolve independently. Thus, to 
determine the evolution of the gi we have to define these couplings more precisely. 
Since our physical input parameters are the masses of the Z and W bosons, it is 
appropriate to define g\ and #2 through the interaction 

£ = \ [(fl?) 2 + (# 2 ) 2 ] (g^W+W^ + \G^V^f) (3.6) 

where V\ = (W 3 , B) are the neutral SU(2)l and U(1)y gauge fields. The bound- 
ary conditions for the coupling constants G ± and G 1 ^ above the SU(2)l x U(1) y 
breaking scale m t are 

G ± = gl, G lJ = Mj with & = (gi, -52) • (3.7) 

After SU(2)l x U(1) y symmetry breaking, the set of /9-functions becomes much 
larger and more complicated. However, if we only work to first order in perturba- 
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tion theory below mt, it is not necessary to derive a full set of /3-functions. The 
only ones needed, (3q± and /3qh, are given in Appendix B. The gauge boson masses 
are related to the coupling constants by 

mw = \v 2 G ± {m z ) , m| = \v 2 tr G^(m z ) . (3.8) 

One can solve the RGEs of Appendix B for G^ with the ansatz G y (t) = g~i(t)g~j(t) 
and therefore G lJ (t) is a rank one matrix for arbitrary t. Thus the neutral gauge 
boson mass matrix maintains its zero mass eigenvalue corresponding to the massless 
photon even below the top quark-threshold. Therefore, it is convenient to define 
gi and #2 in terms of G ± (mz) and G % i(mw) 



We now have to choose a value for MgusY- Of course the simplest possi- 
bility is one in which all supersymmetric particle masses are of order MgusY- 
For completeness we briefly discuss the case of multiple mass scales in the super- 
symmetric particle sector. In the higgsino/gaugino sector we have two free mass 
parameters fi and Mi (it is common to fix M\ by the grand unification relation 
Mi = | tan 2 dy/Mi). We have computed the contributions of the gauginos and the 
higgsinos to P\ t (i = 1, ..,4) and f3 g 2 (i = 1,2). It is clear that the contributions 

of W, B and H to the (3 g 2 can be computed separately and have to be included 

at a scale y/s > M2, Mi and respectively. In the case of the (3^, both gaug- 
inos and higgsinos have to be present at the scale y/s to yield a contribution. As 
a result, the supersymmetric contributions to f3\ i that are proportional to g%, g\ 
or the product g\gi have to be included if the scale y/s > fi2 = max{|/x| , M2}, 
y/s > /i\ = max{|yu| , Mi} or y/s > //12 = maxj/Ji,^}, respectively. The gluino 
contributes only to f3 g 2 at one-loop. The gluino mass (M3) is typically fixed by 
the grand unification relation M3 = (g^/ ' g^jM^- Squark mass parameters are de- 
fined explicitly in section 6. For simplicity we neglect the possibility of generation 
mixing in the squark mass matrices. The result for and f3 g 2 for the MSSM are 
presented in Appendix A. 

We can now compute the effective low-energy coupling constants \i(mz) of 
the Higgs potential by evolving the coupling constants in eq. (3.1) from MgusY to 



I That is, we do not solve the full set of RGEs below mt- For numerical purposes, it is 

certainly sufficient to isolate the one-loop leading ln(m|/m^) terms. 
§ We choose a convention in which M 2 is positive. 
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mz by means of eq. (3.4). Plugging those Aj into eq. (2.10) we find the elements 
of the CP-even Higgs mass matrix. We focus our attention on the simplest case 
where all supersymmetric mass parameters are roughly degenerate. With the first 
order leading log results for the Aj presented in Appendix C these matrix elements 
become in this case 



M 2 n = 



2 2, 2 2 i 

m A s p + m z c p + 



gWz c l 

96^c 2 w 



+ (l2N c 



ml 



-QN C 



mi 



Ml,= 



,2 „2 



m\4 + mis} + [(p f + P g + P 2H ^ In 




+ 12iV c 



m u 
l Z S (3 



Mj 2 = - 



S P C A m 2 A + m 2 z + 



m u 

2 2 
92 m Z 

96^c 2 w 



miySa I \ mi J 



+ - 37V C 



m b 
m 2 7 c 2 



Z^p 




M 2 \ 
iM SUSY ^ 



ml 




(3.10) 



where the constants P,- are 



NJl 



Pf = 



N g {N c [2 - 4s 2 w + 8(e 2 d + e 2 u )s A w ] + [2 - As 2 w + 8s 4 w ]} - P t , 



-44 + 106s 



w 



62s 



w > 



P' = 

9 ~ 

P2H = 

p> = 

r 2H — 



10 + Us 2 w - 26s^ , 
-10 + 2s^-2s^, 
8 - 22s^ + 10s^ . 



(3.11) 



Here the subscripts t, /, g and 2H correspond to the contributions from the top- 
quark, the fermions (excluding the top-quark), the gauge bosons and the two Higgs 
doublets. By diagonalizing the CP-even mass matrix given by eq. (3.10) one obtains 
the individual neutral CP-even Higgs masses and mixing angle a [eq. (2.12)]. With 
the angle a in hand, one then obtains the couplings of the Higgs bosons to gauge 
bosons which are proportional to either sin(/3— a) or cos(/3— a). Finally, the \i(mz) 
and the angle a determine the Higgs self couplings [eq. (2.13)]. These results were 
previously presented in ref. [13]. Similar results have also recently been obtained in 
ref. [17] using an effective potential computation. However, our results differ from 



14 



those in ref. [17] by terms of order g\m\sy V \Yi.{m^ / 'm 2 z ) due to our more precise 
treatment of the theory below the top-threshold [as explained Appendix B] . 

Of course the CP-even Higgs mass matrix may be computed numerically by 
employing the \i(mz) obtained through numerical solution of the RGEs. The re- 
sulting Higgs mass matrix will then be the RGE-improved version of eq. (3.10), in- 
corporating leading logarithmic effects beyond one-loop order. In section 8 we will 
compare these two results to see the numerical implications of RGE-improvement. 

Radiative corrections to the charged Higgs mass sum rule have been obtained 
in refs. [38-40]. For completeness, we also give the one-loop leading logarithmic 
expression for the charged Higgs mass. From eq. (2.9) and (C.6), we obtain 



2 2,2 N c g 2 

m H ± = m A + m w + 



32tt 2 to^ 



2m t m l m 2 (m 2 t m b * 2 4 

—2-2 % I -2" + ~2 ) + 3 m W 

b l3 c p V b l3 c p 



+ i N c(N g -1) + N g + \N H - 10 + 15 tan 2 9 W ] In (^F) , 

W (3-12) 

where the number of Higgs doublets is Njj = 2. The coefficient of l^M^^y/m 2 ) 
in eq. (3.12) is consistent with the exact one-loop calculations of refs. [39] and 
[40]. Note that this result differs slightly from the one quoted in ref. [12], which is 
based on an effective potential computation. This is not surprising in light of the 
remarks made in the Introduction. As discussed in refs. [15] and [39], the effective 
potential is determined from Green functions evaluated at zero external momenta, 
while physical masses are determined by the pole of the (radiatively corrected) 
propagator. This implies that for a particle whose tree-level mass is nonzero, the 
effective potential does not yield the complete leading logarithmic contribution to 
the particle mass. 

The one-loop leading log results [eqs. (3.10)-(3.12)] are useful approximations 
(in the absence of large squark mixing — see section 6) to the fully integrated RGE 
results. This will be discussed more fully in section 8. One clarification is necessary. 
Consider the fully integrated RGE result for some parameter pj. Technically, we 
have only determined pi down to y^i = mt- For ^/s < mt, we have only used 
the one-loop approximation [as in eq. (3.5)] to evolve Pi(s) all the way down to 
\/s = M wea k. Since mt cannot be very much larger than mz, this procedure is 
certainly sufficient for our purposes. 
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4. The Large ra^o Limit 



Until now we have assumed that m^o ~ 0(mz)- In particular, the low-energy 
effective theory (at the electroweak scale) contains two Higgs doublets. Conse- 
quently, when we refer to the VEVs V{ in section 3, we mean Viimz)- In this 
section we shall generalize our analysis to arbitrary m^o. First we consider what 
happens when m^o ^> mz- If we expand eq. (2.11) in powers of m 2 z /m 2 A0 we find 



m 



h° 



= clM 2 n + s\M\ 2 + 2spcpM{ 2 + O 



mi 



m 



A" 



2 2 



9 2 
92 m Z 



w 



mi 



- 6N c c 2 f3-t + 4pP f 



m 



+ (P g + Pm){4 + 4) - 2s 2 c 2 (P' g + P' 2H ) 



In 



^SUSY A 



+ 



4 2 
777. 777- 

12^^ + 6^02/3^- 

m% 



In 



m, 



(4.1) 



+ o 



m, 



m 



.4" 



+ o 



4 2 i / ^Fl A 

^ 2 m z In 



7717 



after using the results of eqs. (3.10) and (3.11). Since eqs. (3.10) and (3.11) 
were derived under the assumption that m^o & 0(mz) it follows that the terms 
proportional to P 2 h and P 2H will be modified when m^o ^> mz by terms of 
O[g 2 \n\\n(m 2 A0 / m 2 z )\ as indicated above. 

Alternatively, we may investigate the case of large m^o by integrating out the 
heavy Higgs doublet. In this scenario one of the mass eigenvalues of A4 2 j is much 
larger than the weak scale. Then, in order to obtain the effective Lagrangian 
at Mweak, we first have to run the various coupling constants to the threshold 
777^40. Then we diagonalize the Higgs mass matrix and express the Lagrangian in 
terms of the mass eigenstates. Notice that in this case the mass eigenstate h° 
is directly related to the field with the non-zero VEV [i.e., f3(m^o) = a(m^o) + 
ir/2 + O(m 2 z /m 2 A0 )}. Below 777^40 there remains only the SM Higgs doublet <p = 
cg$i + s j g$2- The potential is 



V = 7772(0t 0)+ l A(0 t 0) 2 ) (42) 
16 



where the boundary condition for A at m A o is 



\(m A o) = cp\\ + s^A2 + 2s 2 f 3 c 2 fi \z + AcpSpXe + kcpSpM (m A o) 



IWi +9l)c%] (^°) + 384^2.4 ln 



X 



\2N r 



w 



-4 + -J- +6iV c c 2/ 3 -I- 
m~ m% / \m 7 



M 2 



m 



(4.3) 



nv z/ \" L z ^^z 

+ 4 P {p t + p f ) + (4 + 4)(p, + p 2H ) - 2s 2 c 2 (P / + p> 2H) 



where (gf + 5 , 2) c 2/5^ s ^° ^ e evaluated at the scale m A o as indicated. Similarly, the 
mass parameter that appears in eq. (4.2) can be expressed in terms of the soft 
SUSY mass parameters although we will not need this expression to compute m^o. 
The RGE in the SM for A is [41,42] 



16vr 2 /3 A = 6A 2 + | [2g\ + (g 2 + g 2 ) 2 ] -2 £ N Ci h%-\ \^g 2 + \g\ - 2 £ N Ci h 2 fi j , 

(4.4) 

where the summation is over all fermions with hf { = gmf i /(y/2mw)- The RGEs 
for the gauge couplings are obtained from the f3 g 2 given in Appendix A by putting 
Njj = 1. By solving the RGEs for scales mz < \fs < m A o iteratively to first order 
we obtain the light CP-even Higgs mass 



mi = X(m z y = m 2 z c%(m A ») 



+ 



+ 



96tt 2 c 2 



w 



12N, 



c 4 

m% 



6N c c 2f3 



mi 



mi 



+ 4pPt 



In 



rn 



mi 



l2N r 



ml 



6N c c 2 ^ + c 2 w P f + P lg + P 1H 



mr 



mi 



ln 



rn 



A" 



mi 





iwi - 






+ 




+ 6N c c 2 p( 




\ m z 


m%J 



mi 



mi 



3) + 



+ (4 + 4)( p 9 



P 2H ) - 2s 2 c 2 (P' 



P2H) 



ln 



M 2 \ 
iw SUSY \ 



m 



.4" 



o 



m, 



rn 



A°- 



(4.5) 
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Here we have defined 

P 1H = - 94p + (1 - 2s 2 w + 2s 4 w ) c 2 w , 

p ig =4p(- 17 + 70s W ~ 44 ^) - (27 - 36^ + 18s^) 
= ( c /3 + s p)Pg - 2s /3 c /3^ • 



(4.6) 



where the subscripts 1H and lg indicate that these are the Higgs and gauge boson 
contributions in the one-Higgs-doublet model. However, we still must deal with 
implicit scale dependence of c 2 q. Since the fields $j (i = 1,2) change with the 
scale [16,18], it follows that tan/3 scales like the ratio of the two Higgs doublet 
fields, i.e., 



1 dtan 2 /3 d /$ 2 



tan 2 /3 dt ® 2 dt\(Z> 2 



2 =72-71 



(4.7) 



Thus we arrive at the RGE for cos 2(3 in terms of the anomalous dimensions 7, 
given in eq. (A. 10) 



4p( m A°) = cip(m z ) + Ac 2 pc 2 psl(-n - 72) In 



in 



A° 



mi 



(4.8) 



Inserting eq. (4.8) in eq. (4.5), we end up with 



2™2 



2 2 2 / \ , 92 m Z 

™ h o = ™ z c 2p (m z ) + 



w 



ml 



I \,-| 6.Y,r,,^L . ^p f 



+ (P fl + P2ff)(4 + 4) - 2 S 2 C 2 {P' g + P' 2R ) 



+ 



mi 



In 



12N, 



mi 



mi 



m 



mi 



In 



(4 + 4) Pm - ^%s\P' 2H - P m In ( 



^SUSY A 



^SUSY^ 



m 




(4.9) 



which agrees with eq. (4.1). In particular, we have now obtained explicitly the 
term proportional to ln(ra^ ) which accounts for the fact that there are two Higgs 
doublets present at a scale above m^o but only one Higgs doublet below m^o. 
Finally, by taking the limit of large tan (3, we have checked that eq. (4.9) reduces 
to the one- loop leading log result of ref. [5]. 
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5. Physical Definition of tan/? 



In eq. (2.4) we implicitly defined tan/5 in terms of VEVs evaluated at the 
electroweak scale. However, this definition is not appropriate in the case where 
m^o ^> M wca k, since the definition of tan/3 requires the presence of both Higgs 
doublets in our low-energy effective theory. Consider for example the tree-level 
relation for the partial hadronic width of A 

r( A° -> 66) = M^( m 2 -4m2) 1/2 tan 2 /5. (5.1) 

This tree-level relation is true for an arbitrary two-Higgs-doublet model under 
the assumption that down-type (up-type) fermions couple exclusively to Hi (H2). 
Using running parameters evaluated at m^o, eq. (5.1) continues to hold even after 
leading log corrections are included. From eq. (5.1) we can obtain a practical 
definition of tan /3 (771,40). By using the RGE for tan/3 given in eq. (4.7) we can 
obtain tan/3(m^) which at the leading log level matches the definition of tan/5 in 
terms of VEVs given in eq. (2.4). 

It is instructive to show that our results do not depend on the physical definition 
of tan (3. (See ref. [16] for a detailed discussion of the relation between different 
tan/3 definitions.) Consider an alternative definition of tan/5 advocated in ref. [11] 
based on the supersymmetric tree-level relation 

Ami = ml-ml = ml v cos 2(3 + O ( ^) . (5.2) 
L e v \ m zJ 

We can easily obtain the radiative corrections due to the quark/squark sector, for 
arbitrary 771,40. The relevant terms of the Lagrangian are 



Mf + A§ .)] L h + XL ($t L ) ( L t$ .) , (5.3) 



where L = (y,ei) and the indices i,j = 1, 2 run over the two Higgs doublet fields. 
In the MSSM the tree-level quartic coupling constants A and A can be expressed 
in terms of the gauge coupling constants. Clearly, only the terms in eq. (5.3) 
proportional to Kf- cause a selectron-sneutrino mass splitting. The supersymmetric 



* However, there are non-leading log corrections that generate Yukawa couplings of up-type 
(down- type) fermions to Hi (#2)- 
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boundary conditions for h.f- are 

Ag = ^ 2 2 diag(-l,l) . (5.4) 

First, we note that there are no vertex corrections due to the fermions (since 
we can neglect all Yukawa couplings except h\, and ht). Thus the contributions of 
the fermions to the /5-functions are 

f -t f f 167r 2 7/ = -N c h 2 h , , 

vill rema 
heavy Higgs doublet via 



Note that Kf- will remain diagonal at scales below MgusY- After removing the 



$i^ C/ #, $2^S/30, (5.6) 

we obtain Am- = hv 2 A, where A is the coefficient of |-£^0| 2 in the effective scalar 
potential at a scale y/s < m^o. The /3-function and the boundary condition for A 
at the scale m A o are 

A(m A o) = [cgAfi + s 2 M 2 ] (m A o) , 

/ x ( 5 - 7 ) 

p K = -k r f = -A (4 7 f + s}4 ) . 



Thus the resulting slepton squared mass difference is obtained immediately 

Am| = ±v 2 A(M weak ) = m 2 w 



#2 V iW weak / 

V m 40 / V iW wea 



(5.8) 



This equation, which represents the one-loop leading log radiative corrections to 
the slepton mass difference due to the quark/squark sector, is sufficient for our 
purposes. Note that the physical quantities in eq. (5.8) [i.e., mw and Am 2 ] cannot 
depend on the arbitrary scale m A o. It is a simple exercise to check that by taking 
the derivative of eq. (5.8) with respect to \n(m 2 A0 ) we recover the RGE for cos 2/3 
[eq. (4.7)]. 
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6. Radiative Corrections due to Soft Squark Interactions 



In the last two sections we have studied the leading log corrections to the 
Higgs self coupling constants Aj. By inspecting the Lagrangian we find that all 
dimension-four operators of the MSSM respect two U(l) symmetries. Under these 
symmetries, the Higgs fields transform as — > e Mi $j, % — 1,2. One combination 
of these symmetries is gauged [namely U(l)y] whereas the other one is a global 
symmetry [which we call U(l) g ] which imposes constraints on the parameters of 
the theory. One result of the U(l) g symmetry is that the A; (i = 5, 6, 7) remain 
zero even after one- loop leading log corrections are included. However, the U(l) g 
symmetry is broken by dimension-two and dimension-three terms of the MSSM. 
The dominant corrections derive from the squark mixing effects in the top and 
bottom sector. These effects would lead for example to finite (non-logarithmic) 
renormalizations of A5,A6 and A7. In this section, we show how to obtain such 
corrections (see also ref. [8]). 

The most general scalar potential (including Higgs fields and one generation of 
squarks ) takes the following form: 



V° = Vm + Vt + Va + Vq, 



(6.1) 



where we have defined 



V M = (-l) i+j ml$l$j + M| (gtg) + M%U*U + M 2 ~D*D , 



Vr = TP[*lQ)D + TV li$ia 2 Q)U 



u 



■■jr. 



(6.2) 



Aey ( i$i a 2 Qj ) D*U + h.c. 



and Vq contains quartic squark interaction terms. Henceforth, we discard Vq since 
the contributions of these terms to the finite renormalizations of the \ enter only 
at two-loop order. In eq. (6.2), k, I — 1, 2 run over the two Higgs doublet fields, 
the third generation squark fields are defined below eq. (A. 8), and the various Fs 



* Contributions from the sleptons and the other two squark generations are omitted in the 
formulae presented in this section. 
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and A's are determined by the tree-level SUSY relations 

A« = diag [\{g 2 - g 2 Y Q ), h 2 v - \{g 2 - g 2 Y Q )] , 
A Q = diag (h 2 D - \gl \g\ - hfj) , 



A u = diag {-\g\Yu, h 2 v + \g 2 Yu) , 



(6.3) 



A = 



k u = dmg{h 2 D -\g 2 Y D ,\glY D ) , 



(6.4) 



where hjj, hjj, and the squark hypercharges are given in eq. (A. 9) and subsequent 
text, and 

T u = h u (- f i,A u ) , 

r D = h D (A D ,-n) , 

which define the A-parameters: Ajj and Ajj. (The parameter \i also appears in the 
chargino/neutralino sector [43].) The Af k can be expressed easily in terms of the 
Xi (i = 1, 2, ...) of eq. (2.2). If we introduce the collection of fields * = (Q, U*, D*) 
and use the fact that 

d 2 V° d 2 V° 



we obtain the squark mass matrix 



= 



M 2 ab 



9 2 V C 



(6.5) 



(6.6) 



The RGEs due to the squark sector presented in Appendix A can now be 
computed by demanding that the one-loop renormalized potential (in the Landau 
gauge using dimensional reduction [44] and the MS'-scheme) 



V 



32tt 2 



In 



(6.7) 



is independent of the arbitrary renormalization scale a. Here, N c = 3 colors, and 
we have included a factor of 2 since the squark fields are complex. The RGEs for 
the quartic Higgs couplings are 



ik 



dt 



2AgAg + AgAg + AgAg 



\Q\Q 



(6i 



+ AgAg + AgA^ + A (SijStk - S u S jk ) 



The /3-functions at scales below the mass of one or more squark fields are 
obtained from eq. (6.8) by removing the contributions corresponding to these fields 
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and assuming that the coupling constants are continuous. It is clear that the 
dimensionful coupling constants in eq. (6.4) cannot contribute to the /5-functions 
of dimensionless couplings at scales larger than all the mass parameters in V°. 
However, in the presence of dimension-three terms the decoupling of heavy squarks 
becomes non-trivial. To understand what is happening, consider the path integral 
derivation of eq. (6.7), in the case of Mg ^> M^,Mq (i.e., we integrate out D; 
other cases are completely analogous). The generating functional is 



W oc 



dD*dD 



exp <idx 



D*(-i/\)- l D + rf ($*Q)L> + h.c. + C 



(6-9) 

where all other terms of the Lagrangian are included in The inverse propagator 
in the presence of non-zero Higgs fields is 



HA)- 1 = □ - [M 2 ~ + h 2 D <S>j + \Y Q gl{$l - $!)] . 



(6.10) 



The path integral over D and D* is straightforward with the result 

W ocexp^-i J & A xT? ($*Q){-iA) [rf ($*Q)] * + . (6.11) 



When external momenta are much smaller than Mg, the interaction term in 
eq. (6.11) becomes local and can be absorbed into the scalar potential of the low- 
energy effective theory by redefining 



A?' 



h9. 



Ml 

D 



V D V D 
1 i 1 3 



(6.12) 



These are inserted into the /5-function [eq. (6.8)] which is used to run the Higgs 
self-couplings at scales below M^. The logs arising from such a calculation are 

only of 0[ln(M~/M~)]. However, notice that is no longer diagonal and as a 
result the Xi (i = 5, 6, 7) can become non-zero. This can generate phenomena that 
are absent at the leading log level and may be phenomenologically important in 
some circumstances. 

The dimension-three terms also lead to corrections which have no logarithmic 
dependence on the mass parameters of the model. Consider the case of = 
= = Msusy- In this case the corrections described above are not present. 

Nevertheless, important non- logarithmic corrections to the Af k can arise which 
lead to finite shifts in the Xi (denoted by AA^ below). These can be computed by 
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expanding the effective potential [eq. (6.7)] to fourth order in $ as described in 
Appendix D. There are two types of corrections corresponding to triangle diagrams 
and box diagrams. The results for the triangle diagrams (with two powers of 
trilinear coupling constants Ajj, A D or ji) denoted by a superscript (3) are 

AA ( 3 ) _ N c f 2 2 / 2 g 2 2 + g\ \ 2 2 g\ + g\ 
2 16tt2M s 2 usy \ " 17 4/4 

+ f [K-^ 2 )^ + (4-^ 2 )^]} 



AAf } = 



(6.13) 



6 327T 2 M2 USY \ d d \ 4 c y 4 

AA (3) ^ ,2 ( 92 + 9i 2h2 \ , ,2 9l + 9i 

and the results for the box diagrams (with four powers of A^, A^, or ji) denoted 
with a superscript (4) are 
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A\ (4) - Nc f A 4 h A + ;/ 4 /) 4 1 



AA< 4) = ^ ( A 4 ht + An 



A 4 4) = - 967t2 ^4 usy {^ 2 44 + ^ A l h D + h uhW - AuA D ) 2 } 

AAi 4) = - 9fa2 ^ c 4 {^44 + ^ A l h D - h uhW ~ A V A D ) 2 } (6.14) 



l (4)__ N c (l- 

Z SUSY 



As " 967T2M 4 1 D D ^ U U 



(4) _ N cf i f 2 4 3 4 ] 



Finally, the self energy diagrams yield corrections to the kinetic term of the 
Higgs fields which have to be absorbed by redefining the Higgs fields 



For example, the contributions to the A'- coming from the trilinear scalar interac- 



*i->$i= {5ij-\Kj)®j- (6-15) 

ontribut: 

tions are given by 

- / U 2 ~uArr\ _ / Al, -LlAr,\ 

. (6.16) 



4, - * 



y 9 67 r2M2 USY 



/ /i 2 "MtA , , 2 / ^ "Md 



If we then express the quartic terms of the potential in terms of the new fields $ 
we obtain 
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\{gl + gl)A' 22 , 

-\gl{A' 11+ A ! 22 ) , 



§(<7? + <7|)4 2 , 



(6.17) 



where we have used the (supersymmetric) tree-level results for the A^ that appeared 
on the right hand side of eqs. (6.17) evaluated at the scale Msusy- 

Combining all the results obtained in this section, we conclude that the appro- 
priate boundary conditions for the Aj at Msusy are obtained by setting 



where Aj(SUSY) are the scalar self-couplings of the unbroken supersymmetric the- 
ory given in eq. (3.1) evaluated at Msusy- The low-energy effective Higgs potential 
is given by V in eq. (2.2) where the Aj appearing there are Aj(M wea k) obtained by 
solving the RGEs [given in Appendix A] subject to the boundary condition given 
in eq. (6.18). Using the coupling constants so obtained, one may compute the mass 
eigenvalues, mixing angles and coupling constants in eq. (2.9)-(2.13) numerically. 
The effective potential formalism (see, e.g., refs. [6,12,17, and 18]) correctly repro- 
duces the terms that arise from AA^ and AA 4 - 4 \ However, the effective potential 

method does not pick up the terms arising from AAp , which derive from wave 
function renormalization. 

In the next section we shall present numerical results of this procedure. How- 
ever, it is instructive to first look at a few analytic results. In the one-loop logarith- 
mic approximation, the RGEs are solved analytically as shown in section 3. The 
effect of the new boundary conditions [eq. (6.18)] at the one-loop level is simply 
additive. We shall include only the effects of the third generation squark mixing. If 
we denote the one-loop leading log squared mass shifts obtained from eq. (3.10) by 
(Am 2 ) ill, then we obtain the following expressions for the neutral Higgs masses 



Ai (Msusy) = A* (SUSY) + AA^ + A\f> + AX\ 



(6.18) 



2G 



in the limit of large tan (3: 
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assuming that m A o > mz [if m A o < mz then interchange m^o and m#o]. Taking 
the limit (5 — > n/2 on the right hand side of eqs. (6.19) and (6.20) is subtle because 
of factors of Cp in the denominator. The appropriate limit is one where the Yukawa 
coupling h b = g2m b /(\^2m\ycp) is fixed. Thus, in the limit (3 — > ir/2, it follows 
that m b — > such that m b /cp is fixed. We also refrain from setting = 1 to 
preserve the symmetry of the formulae above. 

As a second example consider the radiative corrections to charged Higgs mass. 
Define 



Am 2 H± 



m 2 H± 



m 2 A o 



m 



W ' 



(6.21) 



so that Am 2 H ± 
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where (Am^ ± )iLL is the value of Am 2 H± obtained from eq. (3.12). This result is 
consistent with the one-loop calculations of refs. [39,40] in the limit of large Msusy- 
[In the same limit, the effective potential computations of refs. [12] and [18] differ 
slightly from the above results for the reasons mentioned below eq. (3.12).] Under 
the assumption that all the soft-supersymmetry-breaking parameters are of the 
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same order and tan/3 <C mt/rrn,, the dominant contribution to Amj I± is 

Aml ± = - Nc ^f ( -JLJ] 2 ( ^L_) 4 + o (gWt) ■ (6-23) 
H 64tt 2 V^susy; \spmwj V2t) 

For sufficiently large /i, this correction dominates the leading log contributions 
which grow only as mf . 

7. Numerical Results 

In this section we evaluate the radiative corrections to the Higgs masses and 
couplings. These have been computed by numerically solving the RGEs for the 
Higgs self-coupling parameters to determine the Aj(M wea k) as described in section 
3. These results are then inserted into eq. (2.10)-(2.13) to obtain the radiatively 
corrected Higgs masses and couplings. In all numerical results presented in this 
paper, we shall take the squark mass parameters to be equal to a common soft- 
supersymmetry breaking mass = = = Msusy- 

In fig. 1(a) and (b) we plot the light CP-even Higgs mass as a function of tan/3 
for mt = 150 and 200 GeV for various choices of m^o. All A-parameters and fx 
are set equal to zero. The Higgs mass saturates at a maximum value, m™ ax , when 
tan/3 and m^o become large. Furthermore, m^o converges to 77140 in the limit 
tan/3 — > oo, as long as 77240 < m™f x . The reason for the h°-A° mass degeneracy 
in this limit is easily understood. The tan/3 — * oo limit can be implemented by 
setting m\ 2 = 0. In this case, the model possesses an unbroken global U(l) g 
symmetry which guarantees that m^o = m^o to all orders in perturbation theory. 
That is, the radiative corrections to m^o in this particular limit vanish exactly. In 
the opposite case where m^o > m^, ax , m#o = 777,40 to all orders in perturbation 
theory. In contrast, the radiative corrections to m\ Q are substantial and grow with 
mf. Moreover, it is easy to see that for 777,40 mz, the dominant m^-contribution 
to m^ is independent of tan/3 [see eq. (4.9)]. As a result, at fixed m^o > m^ ax , 
the radiatively corrected m^o will reach a maximum (minimum) at tan/3 ~ oo 
(tan/3 ~ 1), due to the tree-level behavior of m\ on tan/3. 

For fixed tan/3, m^o reaches its minimum value, m™ n , when 777,40 — > 0. Note 
that in contrast to the tree- level behavior (where m^o < 77140), the Higgs mass does 
not vanish as 777.40 — > 0. Moreover, m™ m increases as tan/3 decreases but exhibits 
only a moderate dependence on m t and MgusY- This behavior can be understood 



* If both 4 ^ and /1 ^ 0, then the U(l) g symmetry is not exact, and non-leading log 
corrections to m h a can be generated. 
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as follows. For m^o <C mz and for values of mt and MgusY sufficiently large 
(say, m t <: 2mzsp and Mstjsy ~ 500 GeV), the CP-even squared mass matrix 
[eq. (2.10)] is dominated by the matrix element M22 due to the m\ dependence of 
A2. This yields 

which is in good agreement with the results of fig. 1. One interesting phenomeno- 
logical consequence is that m^o can be larger than 2m^o. This permits a new 
decay-mode h° — > A A which is kinematically forbidden at tree-level. For a de- 
tailed analysis of the h° — > A A decay-mode and its implications see refs. [12] and 
[23]. 

In the limit m^o — > 00 the couplings of h° to gauge bosons and matter fields 
are identical to the Higgs couplings of the SM so that the Higgs sector of the two 
models cannot be phenomenologically distinguished. However, SUSY does impose 
constraints on the quartic Higgs self-coupling at the scale Mstjsy ; an d this does 
influence the possible values of m^o. To illustrate this point, we plot in fig. 2 the 
range of allowed m^o in the case of large m^o (in our plots we take m^o = 300 
GeV). As we have shown above, the lower limit for m^o is attained if tan/3 ~ 1 
and the upper limit is attained in the limit of large tan (5 (in our graphs we take 
tan (5 = 20)} Suppose the top quark mass is known and that h° is discovered with 
SM couplings. If m h o does not lie in the allowed mass regime displayed in fig. 2, we 
could conclude that the MSSM is ruled out. Note that one can also derive upper 
and lower Higgs mass bounds in the SM (at fixed mt) as a function of A. Here, A 
is some high energy scale, below which all Yukawa and Higgs self-coupling A are 
finite (and A > for stability of the electroweak vacuum). The lower SM Higgs 
mass bound is about the same as the corresponding bound exhibited in fig. 2 for 
A = MgusYj while the upper SM Higgs mass bound can be substantially larger 
than the ones exhibited in fig. 2 (if A is significantly smaller than the Planck scale). 
The reason behind this result is the fact that in the MSSM, A is very small (of 
order g%) at A = M SUSY . 

In fig. 3(a)-(d) we plot the RGE-improved CP-even Higgs masses m^o and m#o 
as functions of m^o for m t = 150 GeV, MgusY — 1 TeV and four choices of tan (5. 
In (a) and (b) m^o exhibits a very weak dependence on m^o as long as tan/3 ^ 1. 
In (d) with large tan (3 = 20 we see the near mass degeneracy of h° (H°) with A 
for m^o ^ mz (m^o ^ mz) while mjjo (m^o) stays constant. This behavior is due 
to the global U(l) g symmetry in the limit m\ 2 — > as pointed out earlier. 



f A second maximum for m h a would arise for very small tan /?; however, this regime is ruled 
out because the top- Yukawa coupling develops a Landau-pole at energy scales below Msusy- 



29 



In fig. 4 we plot contours corresponding to the radiatively corrected values of 
m/jo from 20 to 140 GeV in the mj-m^o plane. Results are presented for Msusy = 1 
TeV and A t = = /i = in the case of tan (3 = 1 and 5. In fig. 5, contours of the 
radiatively corrected m^o are shown in the tan /3-m^o plane. Results are presented 
for M SUSY = 0.5 and 1 TeV and m t = 100, 150 and 200 GeV. 

The Higgs production cross-section in a two-Higgs-doublet model via the pro- 
cess e + e~ — > Z* — > ZH°(Zh°) is suppressed by a factor cos 2 (/3 — a) [sm 2 (/3 — a)] as 
compared to the corresponding cross-sections in the SM. In fig. 6 we plot cos 2 (/3— a) 
as a function of m^o for tan (3 = 0.5, 1, 2 and 20, for A t = A^ = /i 2 = 0, MgusY = 1 
TeV and two choices of m t . The behavior is similar to that of the tree- level result 
in that cos 2 (f3 — a) — > as ra^o becomes large. This is expected since for large 
m^o, all heavy Higgs states decouple, while the hPZZ coupling [which is propor- 
tional to sin(/3 — a)] approaches its SM value. Nevertheless, it is interesting to note 
that cos 2 (/3 — a) approaches more slowly as m t increases (i.e., as the radiative 
corrections become more significant). 

Up until now, we have ignored the effects of squark mixing by setting the 
A-parameters and \i equal to zero. We now examine the implications of squark 
mixing on the results obtained up to this point. In order to do this, we must go 
beyond the leading logarithmic approximation and include the effects of non-zero 
A and /j, as explained in section 6. In fig. 7 we plot the Higgs mass m^o as a 
function of tan/3 for mt = 150 GeV and for two choices of m^o. All A-parameters 
are taken to be equal; the four curves shown correspond to fx — A — 0, 1,2 and 
3 TeV, respectively. The behavior of m^o at large values of tan/3 is noteworthy: 
for m^o ^ mz, we see that m^o decreases monotonically with A. In contrast, in 
the case m^o ^ mz and large tan (3, m^o initially increases, reaches a maximum 
at A ?s \/6Msusy, and then falls off rapidly [8]. These behaviors can be obtained 
immediately from eqs. (6.19) and (6.20). 

The results just illustrated have significant implications for Higgs phenomenol- 
ogy at LEP [29]. In fig. 8 we plot the Higgs mass m^o and the factor sin 2 ((3 — a) 
as functions of tan/5 for m t = 150 GeV and MgusY — 1 TeV. We fix the sum 
m^o + m^o = mz in both plots in order to bound m^o while keeping Z — > A°h° 
kinematically forbidden. Fig. 8(a) displays contours of fixed m^o +mj l o = mz- To 
the left (right) of these contours, Z — > A°h° is kinematically allowed (forbidden). 
In fig. 8(b) we see that the decay Z — > Z*h°, with a rate proportional to sin 2 (/5— a), 
can be sufficiently suppressed in the large tan/5 regime to escape detection. On 
the other hand, the rate for Z — > A Q h° is proportional to cos 2 (/5 — a) which is near 
unity for large tan/3 and m^o ^ mz- Thus, Z — > A°h° would be observed in this 
regime unless it is kinematically forbidden. In the absence of the Higgs discovery 
at LEP, we can therefore conclude that the parameter regime to the left of the re- 
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spective curves (for various choices of /i — A) shown in fig. 8(a) are excluded. On 
the other hand, at large tan/3, the parameter regime to the right of the respective 
curves cannot be ruled out based on current LEP data. In particular, for large 
\i = A (and for large tan/3), the true experimental lower limit on m^o [i.e., the 
dotted curve of fig. 8(a)] can be significantly lower than the quoted Higgs mass 
limits of the LEP detector collaborations [26-28]. 

In the search for A at LEP via Z — > A°h°, the phenomenology depends in 
detail on the decay branching ratios of A and h°. The main impact of the one- loop 
radiative corrections is a shift in the Higgs masses such that the new (and typically 
dominant) decay mode h° — > A A is now permitted. We have already noted 
that the tree-level limit m^o < m ^ can be substantially violated when radiative 
corrections are incorporated. In figs. 9 and 10, we depict the region of parameter 
space where m^o > 2m^o. As before, all A-parameters are taken to be equal; we 
exhibit curves corresponding to various choices of A and /i, for = MgusY = 1 
TeV. The parameter space where h° — > A A is kinematically allowed lies to the 
left of [or within] the displayed curves. 

In the analysis presented in this section, predictions for the (radiatively cor- 
rected) Higgs masses are obtained as a function of the soft SUSY breaking pa- 
rameters. However, not all choices of these parameters lead to physically (or phe- 
nomenologically) sensible results. For example, consider the cases of mt = 150 
GeV, M SUSY = 1 TeV and m A o = 0, 40, 100 or 300 GeV. Then the condition 
myfi > (or m/jo > 40 GeV if we take recent LEP limits at face values) rules out 
the parameter region above the solid, dashed or dot-dashed curves in fig. 11. 

Finally we consider the corrections to the charged Higgs mass sum rule. In 
figs. 12-14 we plot the shift in the charged Higgs squared mass due to radiative cor- 
rections, Amjj ± [see eq. (6.21)], as a function of tan/3 for various choices of A and 
/i. We choose mt = 2mw and Mq = MgusY = 1 TeV. Note that in our RGE anal- 
ysis Amjj ± is independent of m^o. For large values of tan/3 (say, m^tan/3 ^ mw) 
the leading log term proportional to m\m\j X^vj/S«c|) [see eq. (3.12)] dominates 
and causes m#± to increase above its tree-level value (except in certain cases where 
the effects of squark mixing become very significant). For more moderate values 
of tan /3, a different term in the leading log radiative corrections, proportional to 
mf/s'p lowers mjj± below its tree- level value. In addition, there is a negative con- 
tribution to Am? H± proportional to mf due to squark-mixing [see eq. (6.23)]. As a 
result, there exists a regime in the SUSY parameter space where Am 2 H ± < —m^ 
[see figs. 12 and 14] which would imply that mjj± < m^o. For example, for 
mt = 2myy, tan/5 = 1 and /i = 2.5MgusY) we find irijj± ~ m^o. Note that the 
radiative corrections in the neutral Higgs sector can impose constraints on the pa- 
rameter space that depend on m^o (see fig. 11). These constraints rule out the 
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area below and to the right of the solid curves in figs. 12 and 14. 



8. One-Loop Leading Log vs. RGE- Improved Results 

In this section we compare the RGE-improved results (equivalent to summing 
leading logs to all orders in perturbation theory) to the one-loop leading log result 
The one-loop leading log result is obtained as described in section 3 by solving the 
RGEs iteratively to first order [see eq. (3.5)]. Such terms correspond precisely to 
a complete one-loop perturbative calculation (with no RGE-improvement) where 
only the leading log terms are retained. However, there is an ambiguity to this 
procedure when it comes to diagonalizing the CP-even mass matrix [eq. (3.10)]. 
If one is performing a full one- loop computation with no RGE-improvement, then 
one should diagonalize this matrix perturbatively and obtain expressions for the 
masses that only involve terms to first order in the leading logs. On the other 
hand, a better approximation is to diagonalize the one-loop leading log CP-even 
mass matrix exactly. The resulting mass and mixing angle will differ from those 
obtained by perturbative diagonalization of the mass matrix by terms that are 
formally of higher order in perturbation theory. However, this result will be a 
better approximation to the full RGE-improved result, as we shall demonstrate 
below. Similar considerations have been touched upon briefly in ref. [16]. 

We introduce the following notation. The one-loop leading log CP-even mass 
matrix [eq. (3.10)] can be diagonalized exactly to obtain the neutral Higgs masses 
and mixing angle. These results will be denoted by the abbreviation ILL. On the 
other hand, if we diagonalize eq. (3.10) perturbatively and only retain the leading 
logs to first order, we will use the abbreviation 1LLP. For example, from eq. (3.10), 
we find 
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where ( m ^o) tree is the tree-level Higgs mass with tan/3 and m^o as the physical 
input parameters* The CP-even mixing angle a and (m/jo) tre e are the tree-level 



* Of course, one has to use a definition of tan (3 that is consistent with tan (3 — jv\ at the 
leading log level. Eq. (5.1) is an example of such a definition. 
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values denned as: 
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Eq. (8.1) plus non-leading log terms is the result one would obtain by performing 
a full one-loop calculation. 

In fig. 2 we see that the ILL results are somewhat larger than the fully inte- 
grated RGE-improved results. This can be understood by examining the dominant 
contributions to f3\ 2 . Here it is important to note that the self-coupling constant 
A2 becomes large at low energies and cannot be neglected. From eq. (A. 2) we 
approximately have 

167r 2 (dA 2 )/dt ~ Q\\ - 6hf + 6h 2 \ 2 . (8.3) 

This equation has an infrared fixed point when the right-hand side vanishes [i.e., 
h 2 |(1 + v / 5)A2 or m h o m l.llmtsp]. Thus the logarithmic growth of the RGE 
result for m^o flattens out for large MgusY as m^o approaches its fixed point value 
[see fig. 2(c),(d)] [45]. 

In fig. 15 we plot m^o as a function of tan/3 for various values of m^o, mt and 
^SUSY- We contrast the various methods for obtaining the radiatively corrected 
value for m^o. The solid curve corresponds to the RGE-improved result obtained 
by numerically solving the renormalization group equations for the \ as described 
in section 3. The dashed curves (ILL) correspond to the radiatively corrected 
m/jo obtained from exactly diagonalizing the one-loop leading log mass matrix 
[eq. (3.10)] and the dotted curve (1LLP) corresponds to the one-loop perturbative 
result given in eq. (8.1). These results provide one of the main motivations for this 
paper. The comparison of the RGE and ILL masses shows only a small difference 
in the predicted mass in the case of Mgugy — 1 TeV. Note that the one-loop 
perturbative formula for m^o [eq. (8.1)] begins to differ substantially from the 
RGE and ILL results for low and moderate values of tan/3. The comparison of 
ILL and 1LLP shows agreement in the large tan/5 limit as expected. Furthermore, 
in the large m^o limit (in fig. 15 we take m^o = 300 GeV) and arbitrary tan/3, 
the difference between the ILL and 1LLP masses is suppressed by a factor of 
0(171^0 / m^o) and is only significant for large m t and MgusY- However, notice the 
extremely large difference between (m^nx and (m/ l o)iLLP m the case of small 
m^o and small tan/5. In particular, in contrast to the RGE and ILL results, 
(m/jo)iLLP is independent of m^o when tan/5 = 1. 



33 



For values of MgusY much above 1 TeV [e.g., MgusY =10 TeV in fig. 15(c) 
and (d)], we begin to see an appreciable deviation between the RGE and the fLL 
results when m^o > mz- In contrast, for m^o <C m 2 z , the RGE and ILL results 
are roughly the same. This is easily understood — in this case, the lightest Higgs 
boson is dominantly Hi and its squared mass is m? — Since there are no 

large mt corrections to M\\, the radiative corrections to M\i are modest and thus 
exhibit very weak dependence on m t and MgusY- As a result, the RGE and ILL 
results are roughly the same. 

To summarize, we have compared three methods for obtaining the radiatively 
corrected m h o. Clearly, the RGE-improved result should be the most complete of 
the three methods examined. The ILL results obtained by exactly diagonalizing 
eq. (3.10) yields results rather close to the RGE-improved results unless MgusY 



Of the three methods for computing the Higgs mass discussed in this sec- 
tion, the fully integrated RGE analysis provides the best approximation. How 
accurate will such an approximation be? To fully answer this question requires 
a detailed investigation of the non-leading logarithmic terms that have been ne- 
glected in our analysis. In this paper, we have identified the most important terms 
of this type only in the limit of large squark mixing. To make further progress 
requires a complete one- loop computation [5,14,15]. The most accurate values for 
the radiatively corrected Higgs mass would then be obtained by identifying the 
complete set of non-leading logarithmic terms (say by subtracting the results of 
eqs. (3.10)-(3.12) from the corresponding exact one-loop calculation) and adding 
these to our fully-integrated RGE results. By comparing our results to those of a 
full one-loop calculation, we conclude that the effect of the non-leading corrections 
on the neutral CP-even Higgs masses is typically no more than a few GeV. Clearly, 
the relative importance of such terms becomes less significant as MgusY increases. 
We can improve our results slightly by identifying the largest of the non-leading 
logarithmic corrections not yet included. In the case of the CP-even Higgs mass- 
squared matrix, when m t ^> mw, the largest of such terms are ones of 0{g\m%). 
From a full one-loop perturbative computation, 



where -M^ GE is the RGE-improved result [i.e., eq. (2.10) with the Aj obtained by 



* For Msusy = 1 TeV the largest discrepancy between these two methods occurs for large 
mt and m^o 3> mz- For example, for tan/3 = 1, m^o = 300 GeV and mt — 200 GeV, we 
find ( 771^0 )rge = 96.8 GeV while (m^o)iLL = 104.4 GeV. In contrast, for the same set of 
parameters we find (m /l o) 1L LP = 111.5 GeV. 



1 TeV.* 
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numerically solving the RGEs). The shift of the light Higgs mass due to these 
non-leading log corrections is of the order of 1 GeV. The largest non-leading log 
corrections to m 2 H ± are given in ref. [40]. With this final improvement, our formulae 
for radiatively corrected Higgs masses are quite accurate in nearly all regions of 
SUSY parameter space for M SU sy k. 400 GeV. 

9. Conclusions 

We have calculated the dominant radiative corrections to the Higgs sector 
parameters in the MSSM. We have obtained analytic formulae for the one-loop 
leading logarithmic corrections as a function of the various supersymmetric pa- 
rameters. Our analysis also includes the most important non-logarithmic effects 
due to squark mixing. In our numerical analysis, we have focused on the case where 
all the supersymmetric partners are approximately mass degenerate. Summation 
of the leading logarithms to all orders in perturbation theory is achieved by solving 
the RGEs numerically. Non-leading logarithmic corrections due to squark mixing 
effects are included by computing the Higgs 4-point functions at the scale MgusY 
and modifying the supersymmetric boundary condition. 

The most significant contribution to the radiative corrections to m^o grows 
with mf and logarithmically with MgusY- A number of the important phenomeno- 
logical implications of the Higgs radiative corrections have already been obtained 
elsewhere (see e.g., refs. [9,12, and 46]). The MSSM cannot be ruled out if LEP- 
200 fails to discover the Higgs boson with a mass m^o ^ mz- In particular, if 
mt ^ 150 GeV, then m^o <; mz for large m^o and tan/3 if MgusY is sufficiently 
above m t . For small values of m^o 30 GeV) and small values of tan/3 1) 
the light Higgs mass can be larger then 2m^o. In this case the decay h° — > A ^ 
is kinematically allowed and provides the dominant decay mode over a significant 
region of the parameter space. The leading log corrections to mjj± only grow with 
m 2 . However, there are non- logarithmic contributions due to squark mixing effects 
that grow as mj. These effects could (in extreme cases) yield m#± ^ m^o. Large 
squark mixing effects can also lower the experimental lower bounds on m^o in the 
large tan (3 region. 

Numerical comparison of the first order leading log corrections obtained in this 
paper with the full one-loop radiative corrections in the limit of large tan /3 [5,14,15] 
shows agreement within a few percent or better for MgusY ~ 400 GeV. The non- 
leading log terms not included in our analysis are thus of the same order as higher 
order leading log terms which are summed by renormalization group improvement 
and both should be included in any full one-loop calculation. The one-loop leading 
log expressions also serve as a useful check of any complete one-loop computation. 
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APPENDIX A 



Renormalization Group Equations 

The /3-functions in the non- super symmetric two-Higgs-doublet model with the 
discrete symmetry $1 — > — $1 and in the MSSM are well known [41,31,47,48]. Here, 
we generalize these results in two respects. First, we take the most general (CP- 
conserving) two-Higgs-doublet potential [i.e., with no discrete symmetry] given in 
eq. (2.2). Second, we explicitly treat supersymmetric particle contributions to the 
/9-functions in the step-approximation. That is, a particular SUSY contribution is 
omitted for scales below the corresponding SUSY particle mass [49]. We then find 

IGttX = {eA 2 + A| + (A 3 + A 4 ) 2 + A| + 12A| + §[2^ + (g 2 + gj) 2 ]}e z 

1 ' * 1 (A.l) 

+ - \hUg\Y Qi + gl) + foi + gfY^} 0^} 

16n 2 p X2 = {eA 2 + A| + (A 3 + A 4 ) 2 + A 2 + 12A| + |[2^ + (^ + g ff}}e z 

+ J2 N ^ {-^ufiufiz + {\glY Di ) 2 e 5 _ + {h 2 Ut + yiYufe^ 

1 ' 1 (A.2) 

+ [h^ + \hlXg\Y Q , - gl) + \{g\ + &q) 

16tt 2 /? A3 = {(Ai + A 2 )(3A 3 + A 4 ) + 1\\ + \\ + A 2 , + 2\\ + 2A? + 8A 6 A 7 

+l[2g$ + (g 2 2-g 2 i) 2 ]}o z 

+ J2 N ^ {-^hlfivfiz + \g\Y Di (h% - \g 2 Y Dz ) 9^ 

i 

- Wu, (h 2 Ut + \g 2 Y Ut ) e Di + hlhlfi^ 

+ [huA - \hUglYQ, + gl) + \h 2 Di {glY Qi - g 2 2 ) + \{g A 2 - g^)] 6^} 

- hl^H + 9hl0w d B e H - Wb B H ~ 167r2A 3(7l + 72) , 

(A.3) 
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16tt 2 /? A4 = [A 4 (Ai + A 2 + 4A 3 + 2A 4 ) + 4A 2 , + 5A| + 5A? + 2A 6 A 7 + \g\g\\ 6 Z 
+ J> Cj {2hlhlfiu£ z - hl t h%6~6 5t - {h\ - yi)(h 2 Di - \g 2 )6~] 

i 

+ 2gl6~e s - 2g 2 1 g 2 2 6 W 6 S 6 S - 167r 2 A 4 ( 7 i + 72) , 



16tt 2 /? A5 = [A 5 (Ai + A 2 + 4A 3 + 6A4) + 5(A| + A?) +2A 6 A 7 ] B z 
- 16tt 2 A 5 (7i + 72) , 

16n 2 p X(i = [A 6 (6A1 + 3A 3 + 4A 4 + 5A 5 ) + A 7 ( 3A 3 + 2A 4 + A 5 )] 9 Z 

- 8tt 2 A 6 (37i + 72 ) , 

16tt 2 /? A7 = [A 7 (6A 2 + 3A 3 + 4A 4 + 5A 5 ) + A 6 ( 3A 3 + 2A 4 + A 5 )] 6 Z 

- 8tt 2 A 7 ( 7 i + 3 72 ) ; 



(A.4) 
(A.5) 

(A.6) 

(A.7) 



where 



/ s \ { 1 for s > My, 
9 x = e(\nA r ) = l - f (A.8) 

V M 2 X ) \o fors<M 2 , K } 

and Mx is the mass of the field X. To first order in M^ eak /M| USY , we can ignore 
mixing in the supersymmetric mass matrices. In this case, the supersymmetric par- 
ticle spectrum consists of the gluino (g), neutralinos (B, W 3 , H®, H®), charginos 
(W ± , H^, L-type squarks and sleptons (Qj,Q 2 ), and R-type squarks and slep- 
tons (Ui,Di). The corresponding supersymmetric particle masses are given by: 
M~ = |M 3 |, = I Mi I, M~ = |M 2 |, Mg = |//|, the L-type squarks [sleptons] 
are degenerate with mass Mq and the R-type squark [slepton] masses are 
and Mp; . The label % runs over three generations and counts both squarks and 
sleptons; the corresponding fermion partners (i.e., quarks and leptons) are denoted 
by Ui and D{. The up- and down-type Yukawa couplings are proportional to the 
corresponding quark [lepton] masses 

hlJi = ~M ~r ' ^ A = ~7! r ■ ^ 9 > 

The color factor N Ci = 3 [1] when the label % refers to (s)quarks [(s)leptons], and the 
hypercharges of the squarks [sleptons] are Yq. = | [—1], Yd. = | [2] and Y\j i = — | 
[there is no Finally, the anomalous dimensions 7 j (j = 1,2) in the Landau 
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gauge are given by 



d i ^2 
7l = -ln $1 



64^ 



^1 + 3^-4^^4^ d z 
~^w 9 H - 2 9i 9 b 9 h\ > 



(A.10) 



72 = -ln$ 2 



64tt 2 



9^ + 3^-4^^^^)^ 



The /3-functions for the U(1)y, SU(2)l and SU(3) C gauge couplings g\ = 
#2 = # and 5-3 = g s are 



48^ 



i 

+ ±N H e t + N 6 e s y )gt, 



48tt 2 # 



.92 



48^ 



(2^ + 0^.) + ±i\fa0 t + N & 6 S + 4A^0~ - 226 




(A.H) 

where Q = 1 [0] for color triplet [singlet] fermions and their supersymmetric part- 
ners. The number of Higgs doublets is Njj = 2, the number of higgsino doublets is 
N^- = 2, the number of fermionic left-handed triplets is Ny, = 1 (the W) and the 
number of fermionic color octets is N5 = 1 (the gluino g). The /3-functions for g\ 
and g2 are valid for scales ^/s > mt where the electroweak gauge theory is unbro- 
ken; for the /3-functions below top quark-threshold see Appendix B. The boundary 
conditions are g^{m 2 z ) = 4tt/(128c^), 4vr/(128s^) and 1.38 for k = 1,2, 3, respec- 
tively. We take = 0.23 in our numerical work. One can explicitly check that 
eq. (3.3) is satisfied for scales larger than the largest SUSY mass parameter. 

It is important to clarify the validity of the above formulae with explicit 6- 
function treatment of thresholds. If all supersymmetric particle masses are roughly 
degenerate and m t ~ mz, then integration of the above formulae will correctly sum 
all leading logarithmic contributions to all orders in perturbation theory. However, 
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if there are non-degenerate thresholds this is no longer the case. For example, the 
treatment of the top quark threshold is discussed in Appendix B. Below this thresh- 
old the electroweak gauge group is broken and thus the effective low-energy theory 
now has many new terms which require the introduction of new couplings. One 
could work out the complete set of RGEs, including those for all couplings (with 
appropriate boundary conditions at ^/s = mt corresponding to the requirement of 
electroweak gauge symmetry). In practice, such a procedure would be overkill; after 
all, mt is not much larger than mz- In Appendix B we show how to correctly ob- 
tain the one-loop leading logarithmic terms proportional to g\\(\(m\ j m?) . These 
are certainly sufficient for our purposes. 

A similar set of remarks hold for the supersymmetric particle thresholds* Let 
Mgygy be the largest supersymmetry-breaking mass parameter in the MSSM. Then 
for y/s < ^susV' we should in principle consider a new low-energy effective theory 
with vertices that are no longer constrained by SUSY. This would introduce many 
new couplings, and the RGEs for these couplings would be required in order to 
sum logarithmic terms such as \n(M? / Mj) to all orders. (Mi and Mj are masses 
of different supersymmetric particles.) Again, this is much more than we need. 
In deriving the RGEs above we have implicitly assumed supersymmetric relations 
among the various couplings of the theory for all MguiSY < ^ < ^SUSY ( w here 
we have assumed that the mass of the lightest supersymmetric particle, Mg^y 
lies above mz)- Integrating the RGEs listed above will then yield terms of order 
g\ ln(M?/Mj), although these terms will not be correctly summed to all orders.^ 
On the other hand, the set of RGEs corresponding to the (non-supersymmetric) 
two-Higgs-doublet model (i.e., for y/s < M^jjg Y ) is sufficient to sum to all orders 
logarithmic terms proportional to l n (^|uSY/^weak)> where MgusY is some aver- 
age supersymmetry breaking mass. Results obtained in this manner are certainly 
accurate enough for our purposes. 

Thus, to complete the set of RGEs required for the computation of the Xi(mz), 
we need RGEs for the Yukawa couplings corresponding to the (non-supersymmetric) 
two-Higgs-doublet model [41,31,47]. These are given by 



* In the numerical work presented in this paper, all supersymmetric masses are taken to be 
degenerate. Thus the considerations below do not affect these results. 

f Such terms are typically not as significant as the terms proportional to g\ ln(mj / 'm z ) , which 
appear multiplied by powers of m t /mz in the formulae for Higgs masses and couplings. 



16tt 2 /^ 



(9 h 2 , l h 2 o 2 9 2 17 n 2\ h 2 

[\h\ + \h\ + h 2 T - Sg 2 3 - \g\ - l 2 g\) h 2 b , 



(A.12) 



40 



All other Yukawa couplings can be neglected. Note that eq. (A. 12) assumes that 
the Higgs-fermion interactions have the same structure as the ones in the MSSM. 
For completeness we present the /^-functions for the mass parameters (although 
these are not required for any of the applications presented in this paper) 



16tt 2 /3 7 



167T 2 /5 m 2 2 = 



[3771^ Al + 77222 (2A3 + ^) — ^ m l2^&] 

[377122A2 + (2A3 + A4) - 6ml 2 h] 
(A3 + 2A4 + 3A5)mf 2 — ^^6 m n — 3A777122 



7i™n 
72^22 



5(71 +I2)m\ 2 . cr 

(A.13) 



APPENDIX B 



Coupling Constant Evolution Below the Top Quark Threshold 

In this appendix we present the /3-functions for the relevant gauge and Higgs 
self-interaction coupling constants below top quark-threshold. Note that by for- 
mally integrating out the top-quark, we include logarithmic contributions to the 
Higgs mass matrix of Olg^m^ ln(m|/m^)]. However, it is important to emphasize 
that such terms are formally smaller or of similar size to some non-leading log 
terms, e.g., terms of 0{g\mf). [The latter can be included easily according to 
eq. (8.4).] Nevertheless, isolating such logarithmic terms serves as a useful check 
of more complete one-loop computations. 

First consider the /5-functions for the Higgs self-couplings. These determine the 
running of the \; the values of Ai(M wea k) enter the calculations of the neutral and 
charged Higgs parameters. In calculations involving the neutral Higgs sector the 
contributions from the top quark can be removed by setting ht = 0. In the charged 
sector we note that some (but not all) of the pieces proportional to h\ arise from 
t-b loops and should also be removed. For the purpose of the calculations presented 
here, the simpler procedure of setting ht — for the neutral Higgs processes and 
ht = hj, = for charged Higgs processes suffices?" 

Next we look at the /5-functions of the gauge coupling constants. First we note 
that below the top quark-threshold the SU(2)l x U(1)y gauge symmetry is broken 
down to U(1) EM . The coupling constants of the photon to charged particles are 



I In a more precise procedure one would have to write down the most general potential for 
four neutral fields and four charged fields invariant under U(1) EM and determine the bottom 
quark contributions to the /3-function for each coupling constant separately. The boundary 
conditions for these coupling constants are obtained by requiring that the potential reduces 
to eq. (2.2) at a scale y/s = m t . 
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constrained by U(1) EM gauge invariance. In contrast, the couplings of the W and 
Z are now independent parameters which thus can evolve differently from m t to 
mz- To determine the evolution of the gauge couplings, we fix as experimental 
inputs the masses of the W and Z bosons mz and mw which arise through the 
interaction 

£ = \ (W&V#V? + G%W+W»-) H*H? , (B.l) 

where V{ = (W 3 , B) are the neutral SU(2)l and U(l)y gauge fields. At the scale 
mt where the SU(2)l x U(1)y gauge symmetry is restored we impose the boundary 
conditions 

G f = Skl91 (B.2) 
G % l x = hl9i9j with gi = (gi, -g 2 ) . 

In general the /5-functions corresponding to the four-point interactions, eq. (B.l), 
have two types of contributions: the box diagrams and the wave function renormal- 
izations. In the following we consider only diagrams involving the top-quark. It is 
clear that the quark loop contributions to the wave function renormalization of the 
Higgs boson and the box diagrams to cancel since the (3 g 2 in eq. (A. 11) do not 
contain pieces proportional to or h%. Below mt we remove the H~tb interaction 
from the theory so that there are no vertex corrections to from the bottom 
quark. However, there are self energy diagrams to both the gauge bosons and the 
Higgs bosons due to the bottom quark. In the regime where tan (3 mw/m^ (i.e., 
hb <^ 92) we can neglect the Higgs boson self energy diagram. Thus the only other 
diagrams one has to consider are the self energy diagrams for the gauge bosons. 
We assert that 

S = G±hi ' (B.3) 
Gg = G»5 kl . 

Then the RGEs of G y and G differ from the RGEs of the gauge coupling only 
by the top-quark contributions 

(no summation over i). The anomalous dimensions 7^ and 7± are defined as 

dln(W+^) dA^V) 

_ dln(^) _ dA^(^) [ ■ ) 

113 ~ dt dt 
where Ayy (V = Vi, W) are proportional to the term of the self-energies of 
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the gauge bosons and A' vv (p 2 ) = dAyy/dp 2 . The contributions involving the 
top- quark are 



7± 



t _ <?2 



167r2 ' (B.6) 



327T 2 

where af = (g- 2 , jSl) and af = (0, 



t _ _ _}__ ILL , _-R_-R^ 
7y — oo„2 l a « a j + a « a i J ' 



APPENDIX C 

One-Loop Leading Logarithmic Higgs Self-Couplings 

In the analysis presented in this paper, the RGEs of Appendix A are solved 
numerically for mt < \/s < Mgugy and iteratively to one-loop order for M wea k < 
t/s < mt using the RGEs described in Appendix B. However, it is instructive to 
solve the RGEs completely iteratively to one-loop order [see eq. (3.5)]. In obtaining 
the results below, we have carefully removed the top-quark from the low-energy 
effective theory at energy scales below m t as explained in section 3. We find 



lr„2 , „2i/ \ , 92 



Ai(m z ) = \[gi + gi](m z ) + 



384tt 2 ( 



, 1 Qs 2 w (1 - 2s 2 w ) tf - 2As 2 w (1 - s 2 w ) tf 2 

■w 



- (42 - 102^ + 60s^) tf - 4 (1 - 2s 2 w + 2s^) tf - 8 (l - 2s 2 w + s^) tf 



m 



Di 



4 4 

m z c p 



+ 12 



m 



2 „4 



9 9 

m z c p 



.a 



+ 6 



4 4 



-6 



m 



Di 



m 2 z c 2 



[cw + s w( e Ui +e Di )]+l + & Di s 2 w + 4e 2 D .siy ) tf 



+^u/wtu\ + (1 - ^uA + Ae 2 v / W ) t% 



(C.l) 
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1 r„2 , „2i/ \ , #2 



A 2 (m z ) = i\g( + gi](m z ) + 



384tt 2 4 



(1 - 2 S ^) t| x - 24s^ (1 - s^) t| 12 



- (42 - 102s 2 w + 60s^) t| 2 -4(l-2s 2 w + 2s^ w ) tf - 8 (l - 2s^ + s^) tf 

+5> 



6^% - 12^% e Vi s 2 w + 4e^4 ) 



+ 6- 



m 



17< 



4 „4 



6 ^2 [ C VK - 4^ + e Di )] + 1 - 4 ef/i ^ + Ae 2 v / W t% 



m 2 z s 2 



+Ae 2 D / w t% + (1 + Ae Di s 2 w + Ae 2 D / w ) t\ 



(C.2) 
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Hmz) = -\[gl + g!](mz) - 384^ | 6s ^ ( X + 2s wO 4 1 + 24 % i 1 ~ s w) 
+ (30 - 42s^ + 124.) tf - 4 (1 - 2s 2 w + 2s A w ) tf - 8 (l - 2s 2 w + s A w ) t 



w 
z 



TT^U o 12 U I 2 1 2 Z) 

-6 2 2 + 45^6^ s^ef/^jj* + 6 2 2 + ^ s W e Di s w e Di t z ' 



m z s p 



m z Cp 



+ -3 



2% l c W ~ swieut + e Di )]+l- Ae Ut s 2 w + Ae 2 v / W ) t% 
m z Sp 



+ -3 



m 



A: 



m 2 z c 2 



[c 2 w + s 2 w (e Ut + e D .)] + 1 + 4e A s^ + Ae 2 D / w ) tf 



(C.3) 

We remind the reader that A3 = A3 + A4 + A5 (in the leading log approximation 
A5 = 0), the number of colors N Ci = 3 [1] and the electric charges of the quarks 
[leptons] are ejj i = 2/3 [0] and eo i = —1/3 [-1]. In addition, we have introduced 
the following notation 



t x 



, M x 
In ' x 



rrir 



l z ' 



Mx > mz , 



i = 1,2 
i = 3 



(C.4) 
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where M x = |Mi|, |M 2 |, \/i\,M^ , M r , , M„ for X = 5, W, H, Qi, U h D { . Further- 

(cgi Ui Ui 

more, we have defined xi> X2, X12 such that 

M~ x = //i = max{|/i|, |Mi|} 

M~ = // 2 = max{|//|,|M 2 |} (C.5) 
M~ 2 = /ii 2 = max{|/i|, |Mi|, |M 2 |} 

If Mx < rnz for some field X, simply set the corresponding = 0. 

The couplings Ai, A 2 , and A3 evaluated at the scale mz appear in the CP- 
even neutral Higgs mass matrix and couplings. For the charged Higgs mass and 
couplings we require A4 evaluated at myy- We find 

\ 4 (m w ) = -\g\{m w ) - j - 6t| + 6 tan 2 w t% + 24 tan 2 9 w t$ 2 

(C.6) 



-44-84 + ^iV c , 



99 9 9 

r> Uj LJj o Uj q j^i I r\ 

4 22 ° 2 2 ~ 2 2 

m W S /3 c /3 m W S p m W C p 



where the definitions of and tj^. are obtained from eq. (C.4) by replacing mz 
with rriw- Finally, in the leading-log approximation, Aj(s) = for i = 5, 6 and 7 
at all scales y/s. 



APPENDIX D 

One-Loop Squark Contributions to the Scalar Potential 

In this appendix we present the derivation of the Higgs four-point functions 
from the one-loop effective potential. The (one generation) squark mass matrix is 
given by 

M 2 = M 2 M + M r + M\ , 

where 

for X = M, T and A and Vx given in eq. (6.2). In the special case = 
= Mq = MgusY |>e., = I-^susy] we can ex P an d the effective potential 
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[eq. (6.7)] 



V = V u + 



^ Ms « SY 



In 



^SUSY A _ 3 



trl 



+ 2A ^SUSY 



In 



M 2 

iM SUSY 
„ 2 



tr (Ml + Mp) 



(D.2) 



+ In tr (Ml + Mp) 2 + — I— tr (Mp + M 2 ) 3 

\ h 1 J Jiw SUSY 

-T^^ tr ( M r + Ml) 4 + 0($ 6 )}. 

If we keep in mind that Mp (M\) contains one (two) power(s) of <£> we find the 
quartic terms of the potential 



quartic - ('I*; &j 



) 



+ 



32tt 2 



In 



M 



M 2 

iW SUSY 



tr(M^) M\- 



y)tr (Mi) 2 



tr (M^) 4 



12M^ USY 



(D.3) 



The traces can now be computed without diagonalization of the mass matrix. It 
is now straightforward to absorb the one-loop corrections into the tree-level terms 
by redefining the Higgs tree-level coupling constants as indicated in eq. (6.13) and 
(6.14). 

In the more general case of unequal diagonal squark mass parameters, the 
computation of V qua rtic is more complicated. Here, one must compute the eigen- 
values of the squark masses (perturbatively in M^ eak /M| US Y) before taking the 
traces. Details of the more general computation and the resulting shifts in the 
Higgs tree- level couplings can be found in ref. [24]. 
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FIGURE CAPTIONS 



1) RGE-improved Higgs mass m^o as a function of tan/5 for (a) m t = 150 GeV 
and (b) mt = 200 GeV. Various curves correspond to m^o = 0, 20, 50, 100 
and 300 GeV as labeled in the figure. All A-parameters and /i are set equal 
to zero. The light CP-even Higgs mass varies very weakly with m A o for 
m A o > 300 GeV. 

2) The range of allowed Higgs masses in the large m A o limit (in our plots we 
take itia° = 300 GeV). The lower limit corresponds to tan/5 = 1. The upper 
limit corresponds to the limit of large tan (3 (we take tan (3 = 20). In Fig. 2(a) 
and (b) we vary m t and keep Msusy fixed at 1 and 0.5 TeV, respectively. 
In Fig. 2(c) and (d) we vary Msusy an d keep mt fixed at 150 and 200 
GeV, respectively. The solid (dashed) curves in (c) and (d) correspond to 
the computation in which the RGEs are solved numerically (iteratively to 
one-loop order). 

3) RGE-improved Higgs masses m^o and m#o as a function of m A o, for Msusy = 
1 TeV, m t = 150 GeV and for various choices of tan (3. 

4) Contours corresponding to (the radiatively corrected) Higgs mass m h o = 20, 
40, 60, 80, 100, 120 and 140 GeV as a function of mt and m A o. Results are 
presented for Msusy = 1 TeV and At = A^ = fi = 0. We exhibit the cases 
of (a) tan/5 = 1 and (b) tan/5 = 5. 

5) Contours for constant m^o the tan /5-m^o plane in the case At — A^ — /j, — 0. 
Two adjacent graphs corresponding to Msusy =0.5 and 1 TeV, and three 
choices of mt are displayed in both cases. The contour labels are exhibited 
explicitly in the Msusy = 1 TeV graphs. 

6) The factor cos 2 (/5 — a) as a function of m A o for tan/5 = 0.5,1,2 and 20 
(dotted, dashed, dot-dashed and solid curves, respectively). Results are pre- 
sented for Msusy = 1 TeV. We consider the case of (a) m t = 150 GeV and 
(b) m t = 200 GeV. 

7) The neutral Higgs mass m^o as a function of tan/5 for (a) m A o = 50 GeV 
and (b) m A o = 300 GeV, for m t = 150 GeV and = M SU sy = 1 TeV. All 
A-parameters are taken to be equal (A — At — A^) ; the four contours shown 
correspond to fi = A = 0, 1, 2 and 3 TeV, respectively. 

8) As a function of tan/5, we plot (a) the neutral Higgs mass m^o and (b) 
the factor sin 2 (/5 - a) for m t = 150 GeV and = M SUS y = 1 TeV. 
All A-parameters are taken to be equal (A — At — A^); the four curves 
shown correspond to ji — A — 0, 1, 2 and 3 TeV, respectively. The sum 
m A o + to/jo = mz is kept fixed in both plots. 
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9) The range of parameters in the tan /3-mt plane where m^o = 2m ao = 40 
GeV for various choices of /J, and the A-parameters, and Mgugy = = 1 
TeV. 

10) The range of parameters in the tan/3-m^o plane where m^o = 2m a® and 
m t = 150 GeV for various choices of /i and the A-parameters, and MguSY = 
Mq = 1 TeV. 

11) Contours of constant m^o = 0, 40, 100 and 300 GeV in the //-tan/3 plane 
for m t = 150 GeV, = M SU sy = 1 TeV. In (a) and (b) m h o = for two 
choices of the A parameters, respectively. In (c) and (d) m^o = 40 GeV. 

12) The shift in the charged Higgs squared mass, Amjj ± , normalized to m^r as 
a function of /i and tan/3 for mt = 2m\y and = MgusY — 1 TeV. The 
A-parameters are set to zero. In our RGE analysis Am^ ± is independent 
of m^p- In addition, we depict in (a) and (b), solid curves corresponding to 
m^o = 0, which depend on m^o as indicated in (a). 

13) The shift in the charged Higgs squared mass, Am^, normalized to m^ 
as a function of (a) the A-parameter and (b) tan/3 for m t = 2m\y and 
Mq = MgusY = 1 TeV. The parameter // is set to zero. 

14) The shift in the charged Higgs squared mass, Amjj ± , normalized to m^ 
as a function of (a) A — At — — /i and (b) tan/3 for m t = 2mw and 
Mq = MgusY = 1 TeV. In our RGE analysis Am 2 H± is independent of m^o. 
In addition, we depict in (a) and (b), solid curves corresponding to m^o = 0, 
which depend on m^o as indicated in (a). 

15) One-loop corrected Higgs mass m^o as a function of tan/3 for mt = 200 
GeV and 150 GeV [(c), (d)] and M SU sy = 1 TeV and 10 TeV. The solid 
curves correspond to the fully integrated RGE-improved Higgs mass. The 
dashed curves correspond to the Higgs mass (m/jo)iLL derived by exactly 
diagonalizing the one-loop leading log mass matrix [eq. (3.10)]. The dotted 
curves correspond to (m/,,o)iLLP given in eq. (8.1). In each case, the lower 
(upper) curves correspond to m^o = 20 (300) GeV. 
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